^ 



c^ 



Wilson mass dependence of the overlap topological charge density 

Peter J. Moran^'^, Derek B. Leinweber^, J. B. Zhang^'^ 

^Special Research Center for the Subatomic Structure of Matter (CSSM) and 

Department of Physics, University of Adelaide 5005, Australia 

^ZIMP and Department of Physics, Zhejiang University, Hangzhou, 310027, 

People's Republic of China 

^CSIRO Mathematics, Informatics and Statistics, 

Private Bag 33, Clayton South, VIC 3169, Australia 



O 

t> . Abstract 

o 



The dependence of the overlap Dirac operator on the Wilson-mass regulator parameter is stud- 
ied through calculations of the overlap topological charge densities at a variety of Wilson-mass 

^" ■ values, using a Liischer-Weisz gauge action. In this formulation, the Wilson-mass is used in the 

negative mass region and acts as a regulator governing the scale at which the Dirac operator is 
sensitive to topological aspects of the gauge field. We observe a clear dependence on the value 
of the Wilson-mass and demonstrate how these values can be calibrated against a finite number 
O ■ of stout-link smearing sweeps. The overlap topological charge density is also computed using a 

D \ pre- smeared gauge field for the input kernel. We show how applying the overlap operator leads to 

j^ ^ ■ further filtering of the gauge field. The results suggest that the freedom typically associated with 

smearing algorithms, through the variable number of sweeps, also exists in the overlap operator, 

CS| ■ through the variable Wilson-mass parameter. 

> 

T^ ■ Keywords: overlap, fat-link fermions, stout-link, smearing, topological charge density, 

lO ■ topology, vacuum structure 
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• 

o 

^ ' 1. Introduction 

The topological structure of the QCD vacuum has been the subject of many lattice investi- 
k^ ; gallons over the years. Local patterns in topological charge fluctuations represent a significant 

^^ ■ aspect of this structure. Moreover, important physical phenomena such as a large r]' mass, 6 

C^ \ dependence, and possibly spontaneous chiral symmetry breaking are directly related to vacuum 

fluctuations of the topological charge. By the axial anomaly, matrix elements or correlation func- 
tions involving the topological charge density operator q{x) can be related to relevant quantities 
of hadronic phenomenology. 

Lattice QCD enables non-perturbative studies of the strong interaction from first principles, 
and should prove useful for studying the important topological structure of the vacuum. Unfortu- 
nately, obtaining a lattice discretization for studying topology is not completely straightforward. 
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e.g. naively discretizing the topological charge density generally leads to non-integer values 
for the topological charge. Physical hadronic interactions also observe an approximate chiral 
symmetry that is described by the theory of QCD, where in the massless limit, an exact chi- 
ral symmetry is realized. Unfortunately, naive transcriptions of the continuum theory explicitly 
break chiral symmetry at finite lattice spacing a. 
The Wilson Dirac operator 111], 



Dw = ^lYiu^fi- -raA^-hmj , 



(1) 



contains the irrelevant Wilson term, r A^/2, that explicitly breaks chiral symmetry at 0(a) in 
order to remove fermion doublers. This lattice discretisation is often improved through the intro- 
duction of a clover term |2], however issues with chiral symmetry breaking still exist. 

One technique that has recently been used to successfully reproduce the light hadron spec- 
trum |13|], is to filter the gauge links prior to applying the Dirac operator. These types of fermion 
actions are typically referred to as UV-filtered or fat-link actions. The term "fat-link" comes from 
the smeared, i.e. fat, links that are used to construct the Dirac operator. One can smear either all 
links USSItIII], only the irrelevant terms HlMlilllil], or even just the relevant terms ifish . 
Incorporating at least some amount of UV-filtering has been shown to reduce the eff'ects of chiral 
symmetry breaking ||4l|7l|jjillll[l2lll5|]. Unfortunately, there is no firm prescription for deter- 



mining the correct amount of smearing to apply to the gauge background. One must find a bal- 
ance between speeding up convergence of the Dirac operator, reducing chiral symmetry breaking 
eff'ects, and removing short-distance physics from the gauge field. Of course, when using a fixed 
number of smearing sweeps risw, with a constant smearing parameter a, the smearing procedure 
only introduces irrelevant terms to the action. The fat- link action therefore remains in the same 
universality class of QCD. Nevertheless, this freedom, in the number of smearing sweeps that 
can be applied to the gauge field, can sometimes be regarded as a drawback to fat-link fermion 
actions. 

The difficulties with implementing exact chiral symmetry on the lattice are summarized by 
the well known Nielsen-Ninomiya no-go theorem [16]. The no-go theorem forbids the existence 
of a local lattice Dirac operator, with exact chiral symmetry, and is free of doublers. However, 
in 1982, Ginsparg and Wilson Ul\\ showed that the physical effects of chiral symmetry will be 
preserved if one can find a lattice Dirac operator, D, satisfying the Ginsparg-Wilson relation, 

Djs+ysD^aDRjsD, (2) 

where 7? is a local operator. Luscher later showed llisl] that any Z), which is a solution of ^, obeys 
an exact chiral symmetry. A popular solution to the Ginsparg-Wilson relation is the Neuberger 
Dirac operator [il9i.i2QL 



D 



Dw(-m) 



^ 



l^(-m) Dw(-m) 



(3) 



which satisfies Eq. Q with R = 1/m. Here we consider the standard choice of input kernel, 
Dyv(-m), the Wilson Dirac operator with a negative Wilson-mass term. To produce an acceptable 
Dirac operator m must lie in the range < m < 2. For m < there are no massless fermions, 
while for m > 2 doublers appear |21]. Varying the choice of m within the allowed range results in 
a ffow of /^-eigenvalues, and facilitates a scale-dependent fermionic probe of the gauge field 1201 . 
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Any value of m in the range (0, 2) should yield the same continuum behavior 11221 123'1. However, 
simulations are performed at a finite lattice spacing a, and empirical studies prefer m > 0.9 124J. 
The overlap Dirac operator is extremely useful for studies of QCD vacuum structure because 
it satisfies the Atiyah-Singer index theorem, and will always give an exact integer topological 
charge. However, the value is not always unique and depends on the value of the Wilson-mass 
parameter |119l \2S, |25l |26|]. Studies of the topological susceptibility y = {Q^)/V, have also 



observed this dependence y2S> l22l]- In particular, the study of Ref. |27] found that;^ varied 
with m for small values of /3, but that this dependence decreased as the continuum limit was 
approached. 

In the following, we extend these previous studies to include an analysis of the topological 
charge density q(x), Q = J d^x q(x), as m is varied. In performing an analysis of the topological 
charge density, rather than;^, we have access to a greater amount of information than that which 
is learnt from the susceptiblity. A change in x can be due to a change in the mean-square of 
the topological charge {q^(x)}, or to a more fundamental shift in the long-range structure of 
the vacuum. As such, it is not possible to understand the underlying change in the topological 
structure from a calculation of ;^. 

A calculation of the topological charge density is also a useful probe of the gauge field, due 
to its strong correlation with low-lying modes of the Dirac operator 1281 129|], which strongly 
influence how quarks propagate through the vacuum. Also, while our focus is on the topological 
charge density, all hadronic observables on the lattice are impacted as we are examining the prop- 
erties of a lattice fermion action. In recent years, the available compute resources and algorithm 
enhancements have reached a point where calculations of q(x) using the overlap operator have 
become feasible I29,.,3Q.. .31,1. 

We visualize the topological density as this is currently the most efl'ective way to view the 
extra information. Our analysis will focus on a comparison between the gluonic topological 
charge density that is calculated following the application of a smearing algorithm (see Sect.O. 
Here our decision is motivated by the growing relevance of fat- link fermion actions. By studying 
difl'erent smearings, we are also able to provide a direct quantitative link to the negative-mass 
Wilson renormalization parameter of overlap fermions. We gain useful insights into the simi- 
larities and difl'erences between these smeared actions and the overlap action, and their relative 
efl'ectiveness for studies of topological vacuum structure. A central conclusion of this study, is 
that the "smoothness" of the gauge field , as seen by the overlap operator, depends on the value 
of the Wilson-mass parameter 

2. Simulation details 

Due to the high computational eff'ort involved in a full calculation of the overlap topologi- 
cal charge density, we consider a single slice of representative 16^ x 32 lattice configurations. 
The configurations were generated using a tadpole improved, plaquette plus rectangle (Luscher- 
Weisz [32]) gauge action through the pseudo-heat-bath algorithm, with/5 = 4.80 giving a lattice 
spacing of a = 0.093 fm. 

Five values of the Wilson-mass in the range (1,2) are used to calculate the overlap topological 
charge density, 

ciUx) = -ir(ys(l-^D)). (4) 
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Figure 1: The overlap topological charge density qov{^) calculated with five choices for the Wilson hopping parameter, 
K. Positive regions of topological charge are colored red to yellow, and negative regions are shown as blue to green. 
From left to right, we have k = 0.23, 0.21, and 0.19 on the first row, with 0.18, and 0.17 on the second. There is a clear 
dependence on the value of k used, with larger values revealing a greater amount of topological charge density. 



Results are reported in terms of the input parameter k, which at tree level is related to m by 



1 



K = 



2(-m)a-\- 8r 



(5) 



with the standard choice r = 1. Note that the allowed range for /c is 1/8 < a: < 1/4, and in the 
interacting theory renormalization leads one to consider 1/6 < /c < 1/4. A single calculation of 
qov(x) for one time-slice will contain 16^ x 32 = 8192 sites of information that must be analyzed, 
and this most easily achieved through direct visualizations. In all figures, we represent regions 
of positive topological charge density by the color red fading to yellow, for large to small qov(x) 
respectively. Similarly, regions of negative topological charge are colored blue fading to green. 
A cutoff' is applied to the topological charge density, below which no charge is rendered. This 
allows one to observe the underlying structure of the field. 



3. Dependence on the Wilson-mass parameter 

The topological charge densities, for the five choices of k, are presented in Fig. [T] A clear 
dependence on k is apparent from the figures, with larger values of a: revealing greater amounts of 
topological charge. This is consistent with expectations since as k is increased the Dirac operator 
becomes more sensitive to smaller topological objects. When using smaller values of k these 
objects will not be felt by the Dirac operator. 



The removal of nontrivial topological objects as k is decreased, bears a striking resemblance 
to the well tested cooling IsJllMlllSlHS and smearing S [H, SO, SH El algorithms. 
In these procedures, the links on the lattice are systematically updated such that the gauge field 
is driven towards a more classical state. This results in a removal of topological charge density, 
as the action is decreased. 

The over-improved stout-link smearing algorithm ll43ll is a modification of the original stout- 
link algorithm [42]. Instead of the standard single plaquette, a combination of plaquettes and 
rectangles are used, with the ratio between the two tuned to preserve topology. In every sweep 
through the lattice, all links are replaced by the smeared links tJ^{x) ll42ll 

U^{x) = tx^{iQ^{x))U^{x), (6) 

with 

Q,{x) = k^lix) - Q^(x)) - ^Tr(Qi(x) - ^^(x)) , (7) 



and 



^^(X) : 



y 



Ul(x), (8) 



where 2^y(x) denotes the sum of the plaquette and rectangular staples touching t/^(x) which 
reside in the ^ - v plane. The ratio of plaquette to rectangular staples is controlled by a new 
parameter e |43]. In the following we use the suggested value of 6 = -0.25, which has yielded 
good results in other studies |29, 44] . For the smearing parameter we select a relatively weak 
value of p = 0.01. This should be compared with the maximum value possible for this combina- 
tion of plaquettes and rectangles, p ^ 0.06. Whilst in the standard stout-link smearing algorithm, 
0.1 is the commonly used value. After smearing, the gluonic topological charge density can be 
calculated, 

qsm(x) = ^e,,,^F;^(x) F'/^(x) . (9) 

In order to fairly compare the two definitions for the topological charge density one usually 
applies a multiplicative renormalization to the gluonic qsm(x) II29I1 . 

qsm(x) -^ Zqsm(x). (10) 

This is because after a relatively small amount of smearing the total gluonic topological charge is 
typically non-integer valued due to the presence of quantum field renormalizations. By matching 
to the overlap topological charge density we can alleviate this bias. 

For this study we have a single slice of the topological charge density and thus can not match 
the total topological charge. Instead the renormalization factor is chosen such that the structure 
of the two field densities can be best compared. The best match to the overlap qov(x) is then 
found by calculating, 

min ^ (qov(x) -Zqsm(x)f , (11) 

X 

as the number of smearing sweeps is varied. Two methods for calculating Z are considered; 

• ^calc = Ijx \qov(x)\ I Zx \qsm{x)\ , 

• Zfit, where the renormalization factor is calculated such that (fTTI) is minimized. 
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Table 1: The number of smearing sweeps, Usw, needed to match the overlap topological charge density calculated with 
the listed value of k. The three methods used to find the best match are detailed in the text. 



The first definition is motived by our aim of comparing the structure of the two field densities. 
The second choice was considered to see if the matching could be improved beyond the first defi- 
nition. We also compare with an alternative matching procedure I45j,i4^ in which one calculates, 



^AB 



XaaXbb 



with 



Xab = (1/V) ^ (^aW - qA) (qB(x) - Qb) , 



(12) 



(13) 



where q denotes the mean value of q(x), and in our case qA(x) = qov(x), ^Bi^) = ^smi^)- Here the 
best match is found when Eab is nearest 1 . In this case, the ratio eliminates any dependence on 
the renormalization factor, Z. 

We first consider Zcaic- The overlap topological charge densities, along with the correspond- 
ing best matches, for three choices of k are shown in Fig. [2l We see that as k is decreased, and 
non-trivial topological charge fluctuations are removed, a greater number of smearing sweeps are 
needed in order to recreate the topological charge density. Again this agrees with expectations 
since the overlap operator becomes less sensitive to small objects as k is decreased, and it is these 
objects that are removed by the smearing algorithm. Comparing the diff'erent definitions in Fig.[2l 
shows good agreement in the topological structures revealed. 

The two methods for calculating the renormalization constant Z, together with the values for 
S, are compared in Table [T] As we move down the table there is a monotonically increasing 
trend in the number of sweeps required to match the value of k. We note that despite some 
minor variation in risw, it is possible to correlate the number of sweeps to the value of the Wilson 
hopping parameter. We note that the average renormalization factor Z ~ 0.7, reflecting the fact 
that withp = 0.01 the gauge fields remain rough after ~ 25 sweeps of smearing. The value for H 
remains approximately constant around ~ 0.75, suggesting that after renormalizing the level of 
agreement between the smeared topological charge density and the overlap density is consistent. 



4. UV-filtered overlap 

Let us now consider the eff'ect of evaluating the overlap operator o n a p re- smeared gauge 
field. This is of some relevance to UV-filtered overlap actions |[l5l|47l|48l|49|], in which all links 
of a gauge field are smeared prior to applying the overlap operator. As already seen in Fig. [2l 
applying the overlap operator is in some respects similar to smearing the gauge field. Of interest 
here is whether the overlap operator, acting on a smeared gauge field, will reveal a topological 
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Figure 2: The best smeared matches (right) compared with the overlap topological charge densities (left) in order of 
decreasing k, where qsmi^) is renormalised using Zcaic- Positive regions of topological charge are colored red to yellow, 
and negative regions are shown as blue to green. There is a clear relationship between k and risw, with smaller k values 
requiring a greater number of smearing sweeps to reproduce the topological charge density. 








Figure 3: A comparison of the overlap topological charge density qovi^) computed using k = 0.19 (left), with q^y(x) 
calculated using the same k, on the same configuration, after first applying 25 sweeps of smearing (right). Positive regions 
of topological charge are colored yellow, and negative regions are shown as blue to green. 



charge density close to the input smeared gauge field, or whether further smearing will be needed 
to match the calculated qovi^)- 

To make comparisons clear, we denote the overlap topological charge density, calculated 
using a smeared configuration as input, by q^^(x). We consider the third Wilson-mass, where 
K = 0.19 and the best smeared match was provided by risw = 25. Figure [3] shows the original 
qov(^) along with the new UV-filtered q^^(x). Far less topological charge density is observed in 
the pre-filtered case. Given the previous results, it is clear that a far greater number of smearing 
sweeps will be required to reproduce q(x) using the gluonic definitions. 

Repeating the same calculation as before we find that 45 sweeps of over-improved stout- 
link smearing provides the best match to the overlap topological charge density. A comparison 
between q^^(x) and the smeared qsmi^) is shown in Fig.lH where Zcaic = 0.85. This is approx- 
imately double the original 25 sweeps required to match the overlap topological charge density, 
once again revealing the smoothing aspect of the overlap operator. These results indicate that the 
filtering that occurs in the overlap operator is independent of the input gauge field. 

5. Conclusion 

Using direct visualizations of the topological charge density, we have analyzed the depen- 
dence of the overlap Dirac operator on the Wilson-mass regulator parameter m. As was hinted at 
by previous studies of the topological susceptibility |24, 27], systematic diff'erences appear in the 
topological structure of the gauge field as m is varied. By comparing qovi^) with the gluonic defi- 
nition of the topological charge density, resolved with a topologically stable smearing algorithm, 
a direct correlation between m and the number of sweeps is revealed. Smaller values of k reveal 
topological charge densities that are similar to using a greater number of smearing sweeps. 

From these observations, one can conclude that the "smoothness" of the gauge field, as seen 
by the overlap operator, depends on the value of the Wilson-mass parameter. This is similar 
to fat-link fermion actions in which the smoothness is directly dependent upon the number of 
applied smearing sweeps. These results indicate that the freedom typically associated with fat- 
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Figure 4: The overlap charge density calculated on a configuration filtered by 25 stout-link smearing sweeps, compared 
with qsmi^) after 45 sweeps of smearing. Positive regions of topological charge are colored yellow, and negative regions 
are shown as blue to green. There is a strong correlation between the objects observed. It appears as though the overlap 
operator has again "smoothed" the configuration. 



link fermion actions, through the number of smearing sweeps, is also present in the overlap 
formalism, through the freedom in the Wilson-mass parameter. 

We also considered the application of the overlap operator to a smeared gauge field, which is 
of relevance to UV-filtered overlap actions. We demonstrated that, regardless of the input gauge 
field to the overlap operator, UV-filtering still occurs via the overlap operator. The strength of the 
filtering is of a comparable strength to that of the overlap acting on a hot, unfiltered configuration. 
When creating a UV-filtered overlap action, one must therefore take care to preserve the short- 
distance physics of the gauge field. 

The topological charge density revealed by the overlap operator is similar to that revealed 
after 20 to 30 sweeps of stout-link smearing with smearing parameter p = 0.01, or 2 to 3 sweeps 
at the standard value of p = 0.1. In this light, it is important to continue investigations into 
the extent to which the properties and phenomenology of the overlap operator can be obtained 
through the use of an efficient Wilson-clover action on smeared configurations. 

Future work could also include gauge configurations generated directly using the overlap 
Dirac operator, or possibly with an alternate overlap definition based on staggered fermions |50], 
which may prove more computationally efficient than the usual Wilson-based overlap operator. 
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